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Abstract — A larger set of simultaneous data from the Apollo 12 lunar sur- 
face magnetometer and the Explorer 35 Ames magnetometer are used to con- 
struct a whole-moon hysteresis curve, from which a new value of global lunar 
permeability is determined to he tx - 1.012 ± 0.006. The corresponding glo- 
bal iifduced dipole moment is 2.1 x 10 gauss-cur for typical inducing fields 
of 10"^ gauss in the lunar environment. From the permeability measurement, 


lunar free iron abundance is determined to be 2.5 ± 2.0 wt. ft. Total iron 
abundance (sum of iron in the ferromagnetic and paramagnetic states) is 
calculated for two assumed compositional models of the lunar interior: a 

free iron/orthopyroxene lunar composition and a free iron/olivine composi- 
tion. The overall lunar total iron abundance is determined to be 9*0 ± b. ( 
wt. io. Other lunar models with a small iron core and with a shallow 
iron-rich layer are discussed in light of the measured global permeability. 
Effects on permeability and iron content calculations due to a possible luna 


ionosphere are also considered • 



INTRODUC TI ON 


Theoretical calculations of whole -moon iron abundance have been made 
by several investigators, often based on meteorite compositional models. 
Previous estimates for total iron abundance have generally been lOt by 
weight (Urey, 1962; Reynolds and Summers, 1969; Urey and MacDonald, 1971; 
Wanke et al. , 1973 ). 

Emplacement of a network of magnetometer experiments on the lunar 
surface by Apollo astronauts in 1969-1972 has allowed investigation of 
lunar iron content and magnetic permeability, using in situ magnetic field 
measurements made simultaneously by instruments on the lunar surface and 
in orbit near the moon. Prom these simultaneous measurements whole-moon 
hysteresis curves can be constructed, from which global magnetic permeability 
can be calculated. (Behannon (1968) previously had calculated an upper limit 
of 1.8 for global permeability using measurements of a single lunar orbiting 
magnetometer.) The dual-magnetometer method was first employed by Pyai and 
Parkin (1971), who calculated whole-moon relative magnetic permeability to 
be p, = 1.03 ± 0.13, The uncertainty in this value was later reduced to 
U = 1.029 ± by Parkin et al. (.1973). This measured permeability was 

used to calculate ferromagnetic free iron in the moon to be 5 - U wt .$, 
and total iron in the moon to be 9 ~ ^ wt .$. 

In this paper the earlier work of Parkin et al. (1973) is extended. 

New values of global permeability, free iron, and total iron are determined 
using more data, improved statistical techniques, and better quantitative 
knowledge of the plasma-magnetic field environment of the moon. New hys- 
teresis curves are plotted using a total of 2703 sets of magnetic field 
averages of data measured simultaneously by the Apollo 12 lunar surface 
magnetometer ana the lunar orbiting Explorer 33 magnetometer. Slopes of 
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the hysteresis curves are used to calculate magnetic permeability of the 
moon, from which we calculate free iron abundance in the lunar interior. 

Then we calculate total iron abundance of the moon for assumed free iron/ 

paramagnetic mineral compositions of the lunar interior. Also, we discuss 

> \ 
other lunar models, one with an iron core and another with- an iron-rich 

layer, in light of the measured global lunar permeability. Finally, Russell 

et al. ( 197*0 have recently made permeability calculations using data from 

a single magnetometer, the Apollo 15 subsatellite magnetometer orbiting at 

an altitude about 100 km above the. moon. The results to date indicate the 

possible existence of a lunar ionosphere between the lunar surface and the 

subsatellite orbit. We consider the effects of • such an ionosphere, should 

its existence be proved, upon the lunar magnetic permeability and iron 

abundance results. 

EXPERIMENTAL TECHNIQUE' 

During times when the moon is immersed in an external magnetic field 
which is spatially and temporally uniform, and plasma interaction effects 
are minimized, the total magnetic field E measured at the lunar surface by 
an Apollo magnetometer is expressible as 

B = n H = H + UjtM (l) 

•where H is the external magnetizing field and M is the magnetization field 
induced in the permeable lunar material (see Fig. l). The relative magnetic 
permeability is ^ = 1 + 4jtk-, where k is magnetic susceptibility in emu/cm 0 . 
(In equation (l) the remanent field at the surface site is subtracted out 
for simplicity. ) The total surface field B is measured by a lunar surface 
magnetometer (LSM) located at the Apollo 12 site on the moon (seienographie 



coordinates 23.4° W. longitude, 3.0° S. latitude). LSM instrument proper- 
ties are described in detail by Dyal «t al, (l9?0) and summarized in a 
companion paper in this volume (Dyal et al. , 1974). Simultaneous measure- 
ments of the geomagnetic tail field H are made by. the lunar orbiting 
Explorer 35 Ames magnetometer, which orbits the moon with 0.5 periselene 
and 5 R aposelene at a period of 11.5 hours. Characteristics of the 
Explorer 35 magnetometer are outlined by Sonett et al. ( 19^7 ) • 

To determine iron content and magnetic permeability of the moon, we 
first construct a lunar B-H hysteresis curve using simultaneous measurements 


of the magnetizing field H external to the moon and the total magnetic in- 
duction B on the lunar surface. The slope of the hysteresis curve is then 
determined and the -whole-moon initial permeability is calculated. This 

value is called the "initial" permeability because it is associated with 

-li- 
the B-H hysteresis curve at very small magnetizing field H (order of 10 ' Oe) 

At these small field to lues, the characteristic "S" shape of the hysteresis 

curve degenerates to a straight line through the origin (Ellwood, 1934) . 

Then from the global permeabil.ity, iron abundance in the lunar interior is 

determined as a function of thermal and compositional models of the lunar 

interior . 


MAGNETIC ENVIRONMENT IN THE GEOMAGNETIC TAIL 
In different regions of a lunar orbit (see Fig. 2), the magnetic 
environment of the moon can have distinctly different characteristics. Con- 
ditions desirable for analysis of lunar permeability and iron content exist 
in regions of the geomagnetic tail where the earth’s magnetic field is 
spatially uniform and temporally constant, and effects of plasma currents 
are negligible. 
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‘ The geomagnetic tail is formed due to an interaction of the earth’s 
permanent , .dipole, field with the plasma flowing radially outward from the 
sun at. an average velocity of ~ U 00 km/ sec ; in effect, the earth's field 
is. swept hack into- a cylindrical region (the geomagnetic tail) on the earth’s 
anti solar side. At the distance where the moon’s orbit intersects the tail, 
the field magnitude -is ~ 10 gammas (lO - ^ Oe). Substructure of the tail 
consists of two ’’lobes": the upper or northward lobe has its magnetic field 

pointing roughly toward the earth, whereas the lower lobe field points away 
from the earth. The moon can pass through either or both lobes, depending 
upon, the characteristics of the particular orbit, the geomagnetic dipole 
axis- orientation, and perturbations of the geomagnetic field by solar wind 
pressures. 

Data must be carefully chosen in the geotail so that interaction modes 
and induction modes other than global magnetisation, are minimized. In 
general, it is possible that the total measured surface field B can include 
field. contributions in addition to the geomagnetic field K and the magneti- 
zation field M induced in permeable lunar material. These other possible 
contributing field modes, some of which are very important in regions of 
the lunar orbit outside the geomagnetic tail, are the following: (l) the 

steady .remanent field Bp due to permanently magnetized subsurface materials 
near the Apollo magnetometer site, (2) the toroidal field resulting 'from 
transverse magnetic .(TM) induction, (3) the poloidal field Bp resulting from 
transverse electric (TE) induction, ( 4 ) the field _B associated with dynamic 
interaction between plasma, flow and the above lunar fields, and ( 5 ) the field 
Bp associated with plasma diamagnetic currents. We examine here each of 
these possible, fields to assess its impact on our study of the lunar magneti- 
zation • 



The remanent field B_ will not affect the hysteresis slope measurement. 
Since B p is constant, its components added to equation (l) would, for 
hysteresis curves plotted for the separate field components of B and H, 
simply result in a shift of the hysteresis curve away from the origin with- 
out altering the slope. Indeed, the radial component of.B R at the Apollo 12 
site has been subtracted out of the lunar hysteresis curve presented in 
Fig. 3. 


The toroidal mode would result from currents driven in the lunar- 
interior by a motional electric field V x H, where V is relative velocity 

of the moon with respect to the ambient bulk plasma flow. This type of field 
has not been detected within experimental error by lunar magnetometers (Dyal 
and Parkin, 1971) even in the solar wind, where V ~ 400 km/sec, much faster 
than the moon could move relative to the geomagnetic tail field during qvuet 
times (lunar orbital velocity is 1 km/sec and geomagnetic tail mo „io^al 
velocity is generally ~ JO km/sec (Mxhalov et al. , 1970)). Furthermore, 
is everywhere toroidal to E and tangential to the lunar sun ace (Schubert 
and Schwartz, 1969 ; Sill and Blank, 1970); therefore, contaminating effects 
of this induction mode can be neglected by use of only the radial components 


of B and H in our hysteresis curve analysis. 

The induced poloidal field B results when time -dependent fluctuations 

J? 

of the field external to the moon (3H/3t ) drive eddy currents in the lunar 

interior. B is the dominant induction mode during times when the moon is 

-P 

immersed in the turbulent, rapidly changing solar wind field; however, in 
the limit of low-frequency or small-amplitude driving field fluctuations, 
as in the tail, 3H/3t - 0 and poloidal induction vanishes. Therefore data- 
selection restrictions are placed on peak-to-peak variations of surface and 


external fields 


(to 


ue discussed lawer,' wo el; 


data obviously 
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contaminated by eddy current fields. Furthermore, since global induced 

fields have dipolar symmetry about the direction of 3H/3t, rather than 

H, poloidal contamination would tend to scatter hysteresis curve data. rather 

than change the hysteresis curve slope. 

\ 

When the moon is in the free- streaming solar wind outside the geomag- 
netic tail^the dynamic interaction fields (B^,) caused by the interaction 
between solar wind plasma and the lunar surface fields B_ have been found 
to be ^ 16 gammas at the Apollo 12 site (Dyal et al., 1973)., Plasma stream- 
ing pressure responsible for compression of surface fields is (l/2 )H.h> V 2 , 
where is proton number density, rtip is the proton mass, and V is the 
plasma bulk speed. We assume that the typical plasma speed in the tail is 
characterized by the average magnetopause motions which are about 70 km/sec 
(Mihalov et al., 1970). Deep in the magnetotail lobes the plasma density 
is substantially below the density of O.l/cirr measured in the plasma sheet 

(Eich et al., 1973). Therefore the plasma- field interaction in the geotail 

a 

should be at least 10 times weaker than in the free -streaming solar wind 
and can be neglected. 

Another possible contaminating field is that due to plasma diamagnetism. 
As plasma particles spiral about the magnetic field lines H, their motion 
induces a field (B^) which opposes H. Plasma diamagnetism in the solar 
wind has been measured (Colburn et al. , (1967); Ness et al. , (1967)) as a 
field magnitude change of about 1,5 gammas by comparing Explorer 35 measure- 
ments in the solar wind with those in the plasma void on the antisolar side 
of the moon. Explorer 35 is unable to measure B^ directly in the plasma 
sheet of the geomagnetic tail, however, since there is no well defined 
plasma void created by the noon in the tail comparable to the antisolar 
cavity formed in the solar wind. An Apollo aagAetfittteier on the lunar surface 
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should be shielded (at least partially) from extralunar diamagnetic effects; 
indeed, we have examined extensive data taken when the moon is near the 
neutral sheet, and we have seen differences as much as 2 gammas in fields 
measured by Apollo 12 LSM and Explorer 35 . These field differences are • 
directed such that they oppose the direction of the earth’s field H, charac- 
teristic of- diamagnetic fields. We have found that we can minimise the 
inclusion of plasma diamagnetic fields in our hysteresis curve analysis by 
eliminating all data for which the measured Explorer 35 field magnitude 
|h| < 6>:10~ - 1 Oe. This data selection criterion is applied sc that only data 
from regions deep in the tail lobes are utilized. 

Diamagnetism in the external environment produces the same effect in 
hysteresis-curve analysis as does paramagnetism in the lunar interior, i.e,, 
paramagnetism determined for the moon is only relative to that of the lunar 
environment. Therefore if diamagnetic field data are included in the radial 
B-H hysteresis curve analysis, the measured global permeability will be 
higher than the true lunar permeability. In our earlier work (Parkin et al. , 
1973)5 the calculated permeability was higher than in this paper, due in 
part to inclusion of some plasma sheet data in the analysis. 

MAGNETIC PE PME ABILITY OF THE MOON 

Magnetic permeability and iron abundance of the moon are calculated 
by analysis of magnetization fields induced in the permeable material of 
the moon. When the moon is immersed in an external field it is magnetized; 
the induced magnetization is a function of the distribution of permecble 
material in the interior. Under the assumption that the permeable material 
in the moon is predorainately free iron and iron-bearing minerals, the lunar 
iron abundance can be calculated from the lunar permeability for assumed 
compositional models of the interior. Since the 
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amount of iron present in the lunar interior should be consistent with the 
measured global magnetic permeability, the permeability in effect places a 
constraint on the physical and chemical composition of the moon's interior. , 
In this section we calculate global magnetic permeability of the moon. Lunar 
iron abundance will be determined from the permeability results in the follow- 
ing section. 

Theory 

For the lunar permeability and iron abundance analysis we use a two- 
layer lunar model (see Fig. l). Vie assume that free iron and chemically 
combined iron in the Fe (ferrous) state are responsible for the magnetic 
permeability of the lunar interior, and we model the moon with a homogeneous 
paramagnetic rock matrix (olivine and orthopyroxene models will be used), 
in which free metallic iron is uniformly distributed. Since the suscepti- 
bility of free iron changes several orders of magnitude at the iron Curie 
temperature (T^), a two -layer model is used with ‘the core- shell boundary 

at the Curie isotherm. For R>R,T<T, and for E £ R , T > T . There - 

c c c c 

for 

fore/R > K any free iron will be ferromagnetic, while at greater depths 
where T > the free iron will be paramagnetic. 

For the two- layer lunar permeability model illustrated in Fig. 1 the 
total field at the outer surface of the sphere is expressed (see Parkin 
et al., 1973): 

B = H (1 + 2G) x + H (1-G) y 4- H (1-G) z (2) 

x y z 

where 

(2r,+l)(p^-l) -• X^(r|-l)(2^.j+i) / • 

G ~ 0 V 3 ) 

(2r i +l)(u 1 -’-2) - 2X J (vl)^i-D 
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Here q = p^ and p 2 are relative permeability of the shell and core, 

respectively. Again, are components of the geomagnetic ta^l field 

external to the moon, measured by the lunar orbiting Ames Explorer 35 --ag" 

netometer; B is measured by the Apollo 12 lunar surface magnetometer. The 
> 

permeability exterior to the sphere is p 0 = 1, that of free space; \ = r ' c /\y 

R and R are radius of the core and the moon, respectively. Equation (2) 
c m •• 

expresses the total surface field in a coordinate system wnich -' 1; * ^ origin 
on the lunar surface at an Apollo magnetometer site: x is directed radially 

outward from the lunar surface, and y and z are tangential to the surface, 

directed eastward and northward, respectively. 

A plot of any component of equation (2) will result in a 3 - H hysteresis 
curve. Equation (3) relates the slope of the hysteresis curve to the lunar 
permeability. The average whole-moon permeability p, is calculated from the 
hysteresis-curve slope by setting - Pp - U in equation (o): 

- 1 ' 


G - * 


p + 2 


( 4 ) 


Hysteresis curve data selection criteria 

The most recent data reduction has resulted in the global lunar 
hysteresis curve shown in Fig. 3. This curve has been constructed using a 
total of 2703 simultaneous Apollo 12 -Explorer 35 magnetometer data sets, 
measured during four orbits of the moon through the geomagnetic tail. 

(in the earlier permeability calculation ( Parkin et al., 1973) a smaller 
quantity of 209 data sets was used, with less strict data selection criteria.; 

The present hysteresis curve data have been carefully selected to mini- 
mize the aforementioned contaminating induction and interaction field modes. 
For reasons described in the previous section, contamination from plasma 

het-LSM i* the, extralunar is minimized by eliminauuig 

—D 
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data points for ’which the magnitude of the external magnetizing field 
|H| < 6x10" 5 0e. 

Since poloidal eddy-current induction (Bp) is dependent upon time rate 
of change of the external field (3H/Bt), contamination from the- poloidal 
mode is minimized by restricting variations in the driving field H. Two- 
minute intervals of H and B data (with j H j > 6x10 ^ Oe) are examined using 
a computer program designed to select intervals during which the Explorer 

- " ^ CZ — ^ 

35 and Apollo 12 data peak-to-peak variations are < 2x10 Oe and < 1x10 v Oe, 
respectively. Data which qualify under these criteria are then averaged 
over the two-minute intervals, and the averages are plotted to construct 
hysteresis curves. For likely lunar electrical conductivity profiles (Dyal 
et al. , 1974), the peak-to-peak criteria imposed on data, limit poloidal 
field induction to magnitude <0.7 gamma. We average data over two-minute 
intervals, however, so that the possible error due to poloidal induction 
is much less than 0.7 gamma. Furthermore, since TL, is not a function of 
H but rather of BH/ot, errors due to poloidal induction could add scatter 
to the linear hysteresis curve, but not change its slope. The same data- 
selection process has also been repeated for six-minute intervals of H and 
B data, with similar results. These peak-to-peak criteria have been employed 
both to minimize poloidal field contamination and to insure uniformity of 


the field H over a region larger than the distance between the Apollo 12 
and Explorer 35 for all selected data. 

We note here that similar peak-to-peak criteria could not be applied 
to simultaneous Apollo 15 LSM and Explorer 35 data since by the time Apollo 
15 was deployed if years after the Apollo 12 data were taken, the Explorer 
35 magnetometer data- were contaminated by a spin tone modulation, of period 
~ o.B ran, which caused peak-to-peak oscillations as much as ys. 10 " Oe 
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during quiet tines in the geomagnetic tail* Therefore as yet we have been 
unable to update the Apollo 15 results which we reported earlier (Parkin 
et al., 1973), and we consider only Apollo 12 data in this paper. 

In our calculations of the whole lo.oon permeability we use only results 
from the radial components of B and H data for several reasons. When the 
moon is in the geotail the Apollo 12 LSM location on the lunar surface is 
such that on the average, the geotail field has a radial component at least 
twice as large as either y or z tangential component, i.e., £ 2 H r r/ . 


Also, for the dipolar magnetization field 



2 M hi 


see 


equations (l) ana ( 2 )). Therefore l-i > Iv M , i.e., the resolution of the 
magnetisation field is at least four times greater using the radial com- 


ponent. Furthermore two contaminating field inodes are minimal in the radial 
component* Toroidal fields B ? are not present in the x-conrponent as explained 
in a previous section, and plasma interaction effects (measured m roe solar 
wind) are much lower in the x direction at the Apollo IT site: 


o 


and B /‘B ~ 0.7 (Dyal et al*, 1972). 

Fx / Fy 

Global magnetic permeability and induced dipole moment 

The global lunar hysteresis curve which meets data selection criteria 
discussed in the previous section is shown in Fig. 3* Apollo 12 averages 
of radial (x) components are plotted on the abscissa with simultaneous 
Explcrer 35 averages on the ordinate 5 2703 sets of selected two-minute 
averages are used. The gap in the center of the curve is a result oi the 


data selection criterion which eliminates data measured when the moon is 


near the neutral sheet. The familiar H S M shape of the hysteresis curve 
degenerates at these low field values to a straight line (Fllwood, 193^*) 
intersecting the origin* The x-component of the remanent field at the 
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set a least squares line has been fitted. The best estimate of the slope 
is 1.008 ± 0.004. 

The least-squares result is obtained as follows. Since both Explorer 

35 and Apollo 12 data have errors assumed to be mutually independent and 
> 

Gaussian in nature, neither data, set can be considered to be the independent 
set. For this reason two linear regression analyses have been used, one' 
with B^_ as the independent variable, and the other using as the indepen- 
dent variable. The two regression coefficients have been calculated to be 
1.005 ± 0.002 and 1.011 ± 0.002, respectively; the ±0.002 error limits are 
standard errors. To get the final value we have bisected the angle between 
the two least-squares lines and calculated the slope of the bisector to be 
1.008. The 95$ confidence interval of the regression analyses is ± 0.004; 
therefore our calculated best estimate of the slope is 1.008 ± 0.004. 

Scatter in the Apollo 12 -Explorer 35 data points of Fig. 3 is primarily 
a result of magnetic inhomogeneities, between the moon and Explorer 35 » small 
contributions from eddy current fields, and instrumental noise and offset 
errors in the Apollo and Explorer magnetometers. These error sources may 
introduce small random fluctuations into the data which will not substan- 
tially affect the slope or intercept of the least- squares line. 

From the slope we calculate global relative magnetic permeability of 
the moon to be ]j, = 1.012 ± 0.006 using equations (2) and (4). Eoth extrema 
are greater than 1.0, implying that the moon, as a whole, acts as a para- 
magnetic or weakly ferromagnetic sphere. 

The global induced dipole moment can be calculated from the global 

permeability. The induced dipole moment is expressed as GR^ H, in units of 

8 

gauss-cm^: the hysteresis curve slope = 2G + 1. Using R m = 1.74x10 cm and 


geomagnetic tail field H = J.O Ge, 


the induced dipole moment js determined 
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to be 2.1x10 gauss-cur . 


IRON ABUNDANCE IN THE MOON 


In this section we calculate free iron and total iron abundances in 
the moon -from the global permeability result of the preceding section. 
Calculations are in general dependent upon thermal and compositional models 
of the lunar interior, and are constrained by the known lunar density and 
moment of inertia. 


Theory 

We refer again to the two- layer model of the noon shown m - f ig. 1. 
The core- crust boundary at R = R c is the iron Curie isotherm (assumed to 
be spherically symmetric); its location is a function of the thermal pro- 


file of the lunar interior. In our calculations we assume the moon is 
composed of a homogeneous mineral (orthopyroxene or olivine) of lunar den- 


sity 3.34 gm/em 3 . Free iron grains are dispersed uniformly throu- 
lunar sphere. For R > any free iron is ferromagnetic while at 
depths where T > T , the free iron is paramagnetic. 


hout the 
greater 


In the outer shell there are both ferromagnetic and paramagnetic con- 
tributions to the total magnetic permeability - 1 + The suscept-i- 

bility of the shell is + k^* where k^ is "apparent” ferromagneti 

susceptibility and 1:^ is paramagnetic susceptibility* The j.erroim^ne'D^^ 
component is metallic free iron, assumed to be composed of multidomain, 
noninteracting grains; the paramagnetic component is Fe combined in the 
orthopyroxene or olivine rock matrix* The expected pressures and tempera- 
tures in the outer shell are such that the ferromagnetic susceptibility of 
iron vm.ll not be substantially altered (Bozorth, 1951; Kapitsa , 1955 )• 


The measured ferromagnetic 


susceptibility of the shell material is an 



apparent value which differs from the intrinsic' ferromagnetic susceptibility 

of the iron because of self-demagnetization of the iron grains and the volume 

fraction of iron in the shell. The apparent ferromagnetic susceptibility 

k. is related to the intrinsic susceptibility in _ according to 
ia > -Li 

O.Olp k, 

k. = -s — =nr^ ( 5 ) 

la 1 + N k^ f 

where N is the shape demagnetization factor of the grains ana p is the volume 
percentage of free iron in the lunar material . For spherical iron grains 
N = 4tv/3 theoretically, but experimentally this value is. found to range from 
3 to 4 (ilagata, 1961). We shall use FI = 3.5 in our calculations. 

Unlike the free iron, the paramagnetic iron in the rock matrix has a 
continuous susceptibility across the Curie isotherm. For the combined iron 
the susceptibility varies according to the Langevin relation: 

k ic " nm2 / 3KT 

s 

where K is the Boltzmann constant, T is absolute temperature, n is the number 
of ions per gram, and m is the atomic moment, m is of the order of a few 
Bohr magnetons e ,g . , for the Fe ion m = 5.25 Mp 5.53 M-g (Nagata, 
1961 ). 

For F. < R the lunar material is paramagnetic only, with susceptibility 
c 

k = k + lu : k„ is the contribution of paramagnetic chemically combined 
2 2c. 2a 2c 

iron and kp a is the contribution of free paramagnetic iron above the Curie 
temperature. Again, is only an apparent value and is related to k 2f* 
the paramagnetic susceptibility of free iron, by an equation similar to ( 5/5 
and kp is dependent cn temperature with a relationship analogous to equation 

( 6 ). 
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Compositional anfi therrnal inode la 

Lunar iron abxmdance is determined for t¥0 conipositional models of the 
lunar interior. Recently Eaula et al. (1974) ’nave determined the lunar 
moment of inertia I to be I/MrJ 2 - 0.3952 ± 0.0045, where M and F ffl are . 
lunar mass and radius , respectively* This value is approximately that Oi 
a homogeneous sphere of constant density (l/MRf = 0.400); therefore we assume 
a moon of uniform density p - 3«34 g/cm , i.e., we will use homogeneous 


c ornpo siti onal model s . 

In, the first compositional model we consider the lunar sphere to oe 
composed of orthopyroxene (xFeSiO^ ♦ (l-xjHgSiO^ ), where x is the mole 
fraction of the Fe 2+ phase present. Free iron grains are dispersed uniformly 
throughout the sphere. In a second model we replace the orthopyroxene with 
olivine (xFe^SiO^ . (l-xjMggSiO^ ). Pyroxenes and olivines have been reported 
to be major mineral components of the lunar surface fines and rod' samples 
(Kagata et al.. 1971 ; Kuseman, 1972 ; Weeks, 1972), with combined iron, present 

P i 

as the paramagnetic Fe"‘ i ‘ ion. The ferromagnetic component of lunar samples 
is primarily metallic iron which is sometimes alloyed with small amounts of 
nickel and cobalt (iiag&ta et al. - 1972; Pierce et al., 19/ 1). xhis tree 
iron is thought to be native to the moon (because of its low nickel content/ 


rather than meteoritic in origin (Strangway et al. , 1973)* Orthopyroxene 
and olivine models are consistent with geochemical studies (Urey et al. , 
1971 ; Wood et al. , 1970; Ringwood and Essene, 1970; Green et al. , 1971) and 
geophysical studies (Toksoz, 1974). 

Since the susceptibility of free iron changes several orders of magni- 
tude at the iron Curie temperature (T ), a two- layer model has been used, 
with the core- shell boundary R, at the Curie isotherm (see Figure 10). 
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The Curie isotherm location is determined from the thermal profile used 
for a particular model. Fig. 4 shows the dependence of on depth in the 

moon for hydrostatic equilibrium, superimposed on temperature profiles 

proposed for the moon by several, authors; the increased pressure of the 

' ■■ _•? o 

interior mil decrease the iron Curie point by about 7x10 C per atmosphere 

increase in pressure (Bozorth, 1951)* For the calculations that follow, we 
have constructed three temperature models to span the range of temperature 
profiles in Fig. 4. Presently we use two-layer temperature profiles for 
simplicity. For model profile T x the Curie isotherm is spherically sym- 
metric and located at R /R = 0.9. Shell and core temperatures are 600 °C 
and 1400 °C, respectively. For the model profile the shell is 500 °C, 
and the core. is 1300 °C, while the Curie isotherm boundary is at R c /R m = 0.85. 
Temperatures are 300 °C and 700 °C for shell and core of model profile T^, 

which has R /R = 0.7. 
c m 

Normally the ferromagnetism of free iron is -dependent on pressure as 
well as temperature., At very low magnetizing field strengths such as those 
of the geomagnetic tail, however, the susceptibility of iron is not strongly 
dependent on temperature below the Curie point (3osorth, 1951). Uniaxial 
stress on iron changes its susceptibility (Kern, 196l); however, hydrostatic 
stress should not affect the susceptibility (Kapitsa, 1955) unless, at very 
high pressures, there is a change in volume (Breiner, 19^7 ) • Therefore we 
assume that the susceptibility of uncombined lunar iron is independent of 
pressure. In the outer shell where T < T,, we define intrinsic ferromagnetic 
susceptibility of the free iron to be k^ f - 12 emu/cm^ (Bozorth, 1951), and 
in the core where T > T , free iron intrinsic paramagnetic susceptibility 
is k^ = 2 x IQ' 4 emu/cm^ (Beykowitz and Kneller, 19&9; Bozorth. 1951). 



Global lunar iron abundance 


Using the information described in previous sections we have generated 
the curves shown in Fig. which relate free iron abundance (q) and total 
iron abundance (Q) to hysteresis -curve slope, as follows. The apparent 
ferromagnetic susceptibility of free iron in the shell (l:^) and the apparent 
paramagnetic susceptibility of free iron in the core (k^ a ) are eaG h calcu- 
lated as a function of the free iron abundance (q) using equation (5). The 


shell-core boundary is defined by the. Curie isotherm used for a particular 

pj. 

temperature model. In addition, the mole fraction of the Fe ' pnase in para- 
magnetic rock used in the model ( orthopyroxene or olivine) is constrained 

2 + 

by q and the bulk density. Mole fraction of Fe ' is related to suscepti- 
bility using the experimental data of IJagata et al. (195?) for olivine and 
Akimoto et al. (1958) for orthopyroxene. Furthermore, using the Langevin 

Pj. 

relation (equation (6)), we relate susceptibility of Fe to temperature 
and thus find S's^ ’ susceptibility in the shell ( j and the core 
our three temperature models. Then we combine susceptibilities to obtain 
total shell and core susceptibilities, k n - k^ + k^ c and k^ = + k 2e » 

as a function of q, after which we use k^ and kp in equation (3) to rej.at»e 
free iron abundance q directly to G as shown in Fig. 5* G is related to 
the slope of the hysteresis curve : slope = 2G + 1. 

We complete Fig. 5 by determining Q, the total iron aoundance in the 
moon, as a function of G. Total iron Q is the combined abundance of free 
iron and ferrous iron. From the previously determined mole fraction of 
the Fe^ + phase in the paramagnetic mineral, we calculate the mass of ferrous 
iron as a function of q. Then we add the ferrous iron mass to free iron 
mass and get the total mass of iron in the moon, and thereoy the uotal iron 



We note that for q > 1 wt. % the susceptibility of the moon is dominated 
by the ferromagnetic iron (generally k* » k^ Q >■ k^ ~ k^ a ) and therefore 
the relationship between G and q is independent of the composition, of the 
ferrous component (olivine or orthopyroxene), 'Total iron abundances differ 
for the two compositional mode is because of the different iron contents of 
olivine and orthopyroxene. 

From our measured hysteresis slope of 1.008 1 0.004 we calculate 
G = 0.004 ± 0.002. Using this range of G, we find, from Fig. 5, free iron 
aboundance q and total iron abundance Q, for each thermal and compositional 
model. These iron abundances are summarized in table 1. . . C Table 1 

Fig. 5 shows the ranges of free iron and total iron abundances, which 
are functions of tenperature in the lunar interior and are bounded by errors 
in the hysteresis curve slope. The lunar free iron abundance ranges between 
4.5 and 0.5 wt. these limits correspond to thermal profiles and T v. 
respectively. These thermal profiles have been selected to be upper and 
lower limiting cases, as can be seen upon comparison of Fig. 4 with thermal 
profile descriptions in the text. We accordingly calculate our free iron 


abundance best value as q = 2,5 ±2,0 wt. fo. 

Total iron abundance (Q) in the moon is, in addition, dependent upon 
compositional model. For the free iron/orthopyroxene model, upper and lover 
limits on Q, are 13.7 and 11.8 wt. $, respectively; for free iron/olivine, 
limits are 6.8 and 4.3 wt. ■ Assuming that the moon is composed of one or 
a combination of these minerals and has a total iron abundance between 13.7 
and 4.3 wt. io, we calculate the total iron abundance best value as 
Q = 9,0 ±4.7 wt. Free iron and total iron abundances are shown in 
Fig. 6. 


Fig. 6 
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We note here that the susceptibilities of both olivine and orthopyroxene 


are about an order of magnitude too small to account for the measured lunar 
permeability without some ferromagnetic iron present. Also, we can calcu- 
late the minimum free iron abundance in the moon consistent with the 
hysteresis-curve measurements. To do this, we consider the extreme case 
where the measured whole-moon permeability is assumed to correspond entirely 
to ferromagnetic iron in the outer shell of the moon where the temperature 
is below the Curie point. For this case the bulls lunar iron aomciance 

Q, is 0.9 ± 0.5 wt. c p. This result is independent of Curie-point depth of 

extreme . 

our three models. The/lower limit placed on the lunar free iron abundance 

by our analysis, therefore, is 0.4 wt. io. 


Considerations of an iron core and iron-rich layer 

The whole-moon permeability has also been used to investigate the 
magnetic effects of a hypothetical iron core in the mocn. Density and 
moment of inertia measurements for the moon limit the size of such a core 
to less than 500 tan in radius (Toksoz , 1974). If this hypothetical iron 
core were entirely paramagnetic and the surrounding core were oi ..hCj-.yTOv.Cj.ie 
of average temperature 1100 °C the global permeability would oe 1. 000.5, 

This value is small compared to the measured permeability of 1.012 ± 0.006, 
implying that if such a small paramagnetic iron core exists, its magnetiza- 
tion is masked by magnetic material lying nearer to the surface. Therefore 
the hysteresis measurements can neither confirm nor rule out the existence 


of a small iron core in the moon. 

An iron-rich layer in the moon has been considered by several investi- 
gators (e.g., Wood et al. , 1970; Urey et al,, 1971j Gast and Giuli, 1972^. 
It is possible that early melting and subsequent differention of the outer 
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of a high- density, iron-rich layer beneath a low-density, iron-depleted 
crust. Constraints have been placed on an iron-rich layer by Gast and 
Giuli (1972) using geochemical and geophysical data (for example, measure- 
ments of lunar moments of inertia). One set of their models consists of 
high-density layers between depths of 100 km and 300 km. At a depth of 
100 km the allowed layer thickness is 12 km; the thickness increases with 
increasing depth, to 50 km at 300 km depth. Also presented are a set of 
layers at 500 km depth. By using exactly the same considerations as were 
used in the iron abundance calculations, we calculate whole -moon permeabili- 
ties which would be expected from lunar models with these iron-rich layers. 
The calculations indicate that iron-rich layers allowed by geophysical 
constraints as outlined by Gast and Giuli, if wholly above the iron Curie 
temperature and therefore paramagnetic, would yield global permeabilities 
of about 1.00006. As for the case of a small lunar iron core, the magneti- 
zation field of such paramagnetic layers would be masked by ferromagnetic 
materials elsewhere in the moon, and the hysteresis curve measurement's can 
neither confirm nor rule out these layers. This conclusion would particu- 
larly apply to the Gast-Giuli layers at 500 km depths, which would be almost 
certainly paramagnetic . 

If the iron-rich layers existed shallow enough in the moon to be below 
the Curie temperature and were therefore ferromagnetic, then the measured 
global permeability would be about 3*5* This value is well above the upper 
limit for the actual measured permeability of 1.012 ± 0.006, and therefore 


the Gast-Giuli layers can be ruled out if they are cool enough to be ferro- 
magnetic. It is important to note that the high-density layers discussed 
by Gast and Giuli (1972) can be thought of as limiting cases and that there 

ere innumerable less dense and thinner layers whrcn are allowed 0 j §oopiij?s re 
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geochemical and magnetic constraints. 

E f f e cts of a possible lunar ionosphere on permeability and iron abundance 


calculations 

Measurements of the Apollo.. 15 sub satellite magnetometer in the geomag- 
netic tail have indicated the possibility that an ionosphere exists in the 
region between the lunar surface and the subsatellite 100 tan mean altitude 
(Russell et al. , 1974). Also, results from the Rice University suprathermal 
ion detector experiments indicate that charged particles measured on the luna 
surface could be from a lunar ionosphere (Lindoman et al. , 19/3). If a 
global lunar ionosphere does exist, then it could form a diamagnetic region 
around the moon which could lower the geotail field H to a smaller magne- 
tizing field K* and in turn result in a smaller total surface field B T . 

Then plots of Apollo 12 data versus Explorer 35 ctata would be plots of 
B ’ vs b . causing the measured permeaoility to r.;e lower than ^ru^- 
lunar permeability (in contrast to the plasma diamagnetism eilcct discuses 
earlier, which would tend to male the measured permeability higher). 

If the assumption of Bussell et al. (1974), that a spherically 
symmetric diamagnetic ionosphere fills the entire region between the lunar 
surface and the subsatellite 100 tan altitude, were found to be correct, 
then our global lunar permeability would be adjusted upward, V<e haie 


determined the adjusted permeability and iron abundance values 


by modifying 


equation ( 3 ) for a two- layer (shell/core) permeability model, where in this 
case the shell represents a lunar ionosphere of 100 tan thickness and homo- 


geneous permeability and the core represents a lunar sphere of bulk 
permeability u- T * ne modified version of equation ( 3 ) is 
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( 7 ) 


(2|j <i +p, ) (jjl^+2 ) - 2X^(p, i -^)(|j, i -l) 

where X = R /R : R and R are radius of the core (here the lunar globe) 

' c' s c s 

and the outer radius of the shell (ionosphere), respectively. From the 
results of Russell et al. (1974) we obtain g = - 0.0101 ± 0.0039. Then 
using solutions of equations (3) and (7) we find that our permeability best 
value would be adjusted upward slightly from 1.012 to 1.017, provided a 
lunar ionosphere exists. The corresponding free iron abundance oest value 
would be adjusted upward from 2.5 to 3.9 wt. io. Total iron content best 
values would be adjus'ted downward slightly from 12,8 to 12.1 wt, $'for the 
free iron/ orthopyroxene model, and from 5*5 to 4,7 wt. °jo for the free iron/ 
olivine model. Detailed calculations of the effects of a lunar ionosphere 
on lunar iron content determinations will he deferred until a more detailed 
ionospheric model is presented. 



SUMMARY AND CONCLUSIONS 


(1) Simultaneous measurements by lunar magnetometers on the surface 

of the moon and in orbit around the moon are used to construct a whole-moon 

hysteresis curve, from which the global lunar relative magnetic permeability is 
■> 

determined to be 1 *012 ± 0 .006, 

(2) The corresponding global induced magnetization dipole moment is 

expressed a = 2.1 x 10 22 K. For typical geomagnetic tail fields of H " 10 4 Oe, 

18 3 

the corresponding induced dipole moment is 2.1 x 10 gauss cm . 

(3) Both error limits on the magnetic permeability value are greater than 
1.0, implying that the moon as a whole is paramagnetic and/or weakly ferromag- 
netic. Assuming that the ferromagnetic component is free metallic iron of multi- 
domain, noninteracting grains, the free iron abundance in the moon is calculated 
to be 2.5 ±2.0 wt . %. 

(4) A free iron abundance extreme lower limit of 0.4 wt. - 0 is 

under the assumption that the global susceptibility is due entirely to free iron in 
the ferromagnetic state. This lower limit is independent of composition of the rook 

matrix making up the bulk of the moon. 

(5) Total iron abundance in the moon is determined by combining free 
iron and paramagnetic, iron components for two assumed lunar compositional model 
of orthopyroxene and olivine. For an orthopyroxene moon of overall density 
3.34 g/cm 3 with free iron dispersed uniformly throughout the lunar interior, 

the total iron abundance is 12.8 ± 1.0 wt. %. For a free iron/olivine moon the 
total iron abundance is 5.5 ± 1 .2 wt. %. A summary of iron abundance calcula- 
tions is given in Table 1 . Using extreme upper and lower limits in Table 1, the 
overall total iron, abundance is expressible as 9.0 ± 4.7 wt. %. 

(6) Lunar models with a small iron core and with an iron-rich layer are 

also discussed using the measuroo gjooai lunar permt-aoim^ ■*- — ? , 
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small pure iron core of 500 km radius (the maximum size allowed by lunar 
density and moment of inertia measurements) , which. is hotter than the iron Curie 
point (T>T would. not be resolvable from the data since its magnetization field 
would be small compared to the induced field. we measure, Similarly, an-iroiv 
rich layer in the moon could not be resolved if the iron is paramagnetic, i,e. , the 
iron is above the iron Curie temperature , Gast and-Giuli (1972) have proposed 
a family of high ■'density layer models for the moon which are geochemical J y . 
feasible. If these, models are iron-rich layers lying near the lunar surface so 
that T<T^, the ferromagnetic layers would yield a global permeabiiiiy value well 
above our measured upper limit. Therefore we conclude that such shallow 
iron-rich-layer models are not consistent with our magnetic permeability 
measurements . 
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TABLE 1 


IRON ABUNDANCE OF THE MOON AS A 
OF THERMAL AND COMPOSITIONAL 



Free Iron 

Abundance, wt. fj 

' thermal 
model 

c ornpos i t i o 

model 

T 

3 

i 

i 

T 2 *1 ■ 

orthopyroxene 

1.0 ± 0.5' | 

l i 

i 

2.0 ± 1.0 j 3-0 ± 1.5 
I 

L — _ J 



1 i 

i i 
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o Li. vine 


1.0 ± 0.5 


2.0 ± 1.0 


3.0 ± 1.5 


:tion 
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Total Iron Abundance, wtv % ■ 















FIGURE CAPTIONS 


Fiff. 1 Maenetization induction in the moon. When the moon is immersed in a 

O w 

uniform external field H (in this case the steady geomagnetic tail field) , 

a dipolar magnetization field M_is induced in permeable material, in the 

> 

lunar interior, with the dipole axis of M aligned along the direction of 
H. The total magnetic field near the moon is B = H + 4 n M . The magnetic 
permeabilities of the two layers are p^ and P 2 ' 1 ogions outside 

the moon, p = P Q - 1 (free space) . H is measured by the lunar orbiting 
Explorer 35, whereas B is measured by an Apollo lunar surface magneto- 
meter (LSM) . Measurements of B and H allow construction of a B-.H 
hysteresis curve for the sphere, from which permeaoility and ii on 
abundance can be calculated. 


Fig. 2 


Magnetic environment of the moon during a lunar orbit, with emphasis 


on the geomagnetic tail region. 


The plane of the lunar orbit very nearly 


coincides with the ecliptic plane of the earth's orbit. The earth's 
permanent dipole field is swept back into a cylindrical region known as 
the geomagnetic tail; at the lunar distance the fielo magnitude is 10 
gammas (10~ 4 oersteds) . Also shown is the shock surface ("bow shock") 
created by supersonic flow of solar wind plasma past the geomagnetic 


tail. Substructure of the tail consists of two "lobes' 1 ; the uppei 01 
northward lobe has its magnetic field (dark arrow) pointing roughly 
toward the earth, whereas the lower lobe field points away from the earth 
The moon is immersed in the tail about four days of each orbit; the moo*. 


can pass through either or both lobes, depending upon the characteris 
tics of the particular orbit and the orientation of the earth's magnetic 


dipole axis . 
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Fig. 3 


Fig. 4 


Fig. 5 


Hysteresis curve for the moon. Data points are 2703 simultaneous 
2-minute averages of radial components of the external geomagnetic 
field data H (measured by the lunar orbiting Explorer 35 Ames magneto- 
meter) and total magnetic induction B_ = pH (B is measured by the Apollo 
12 lunar surface magnetometer). Data points are selected from four 
lunations of measurements made when the moon is immersed in the uniform 

geomagnetic tail field, far from the neutral sheet in the tail. In this 

-4 

low-external-field regime (~ 10 gammas or 10 Oe), the hysteresis curve 
is linear and is fitted by a least-squares line of slope 1.C08 ± 0.004. This 
slope corresponds to a whole-moon magnetic permeability of 1.012 ± 0.006. 
The least-squares line intersects the origin exactly in this figure because 
the vertical-axis intercept (the radial or x-component of the remanent 
field at the Apollo 12 site) was subtracted out after the least-squares fit 
was made. 

Temperature profiles for the lunar interior published by various authors: 
(1) Hanks and Anderson (1972); (2) Dyal et al . (1974); (3) Toksdz 
(1974); (4) Sonett et al. (1971). Superimposed is the pressure- 
dependent Curie temperature for iron versus depth in the moon. 

Iron abundance in the moon as a function of global hysteresis curve 
slope. Free iron abundance (a) and total iron, abundance (Q) versus 
the parameter G (hysteresis curve slope equals 2G+1) for three 
temperature profiles described in the text. Total iron abundance is 
shown for two lunar composition models, orthopyroxene and olivine. 
Arrows below the horizontal axis show the range of the parameter G, 
experimentally determined from the hysteresis curve slope: G = 0.004 
± 0.002. The shaded region defines the allowed values of free and total 



iron abundances, bounded by hysteresis curve error limits and by 
thermal models Tj and . 

Fig. 6 Summary of global lunar magnetic permeability and iron abundances 
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